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^ Abstract 
H 

^ We have monitored 8.6 million stars in the Large Magellanic Cloud for 1.1 years and have found 

3 events consistent with gravitational microlensing. We place strong constraints on the Galactic 
halo content in the form of compact lensing objects in the mass range 10~ 4 M to 10~ 4 M . Three 
events is fewer than expected for a standard spherical halo of objects in this mass range, but appears 
to exceed the number expected from known Galactic populations. Fitting a naive spherical halo 
model to our data yields a MACHO fraction / = 0.20±° H, which implies a total MACHO mass 
(inside 50 kpc) of 8.0±g 4 x 10 10 M , and a microlensing optical depth 91^ x 10" 8 (~ 68% CL). 
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There is strong evidence from the observed flat rotation curves of spiral galaxies that such 
galaxies (including our own) have extensive halos of dark matter [1]. We have undertaken an 
experiment that uses gravitational lensing to search for the dark matter in our own Galactic 
halo [2]. The experiment is sensitive to lensing by compact objects over a broad range in mass, 
spanning from ~ 10 -7 to ~ 10 M . This mass regime contains a number of plausible dark matter 
candidates, including brown dwarfs, black holes, and other stellar remnants. Astrophysical objects 
of primordial elemental abundances are thought to require a minimum mass of 0.08 M to ignite 
hydrogen fusion, so this experiment is sensitive to non-luminous objects that would have escaped 
detection in existing sky surveys. Furthermore, we note that a Galactic halo composed entirely of 
baryonic objects is consistent with the nucleosynthesis bounds on cosmic baryon density [3]. 

As suggested by Paczynski [4], a population of MAssive Compact Halo Objects (MACHOs) 
would be detectable by their gravitational 'microlensing' influence on background stars. Microlens- 
ing refers to the special case of gravitational lensing in which the splitting of the source star into 
multiple images is too small to be resolved, but the lensing phenomenon causes a change in the 
apparent brightness of the source which is time-dependent due to the relative motion of the Earth, 
lens and source. The apparent amplification A(t) of a background star, in the point-source point- 
lens approximation, is given by [5] 



1/2 

(1) 

where t = 2te/v is the time for the lens to move through two Einstein radii (r^) relative to the 
line of sight, r E = ^jAGmD/c 2 , m is the lens mass, D = D lens (D SOUIce - D lens ) / D source , and u(t) 
is the distance between the lens and the line of sight in units of te- 

Three teams have reported the detection of gravitational microlensing events, and the world 
total now exceeds 50 events [2,6-9]. Most of these have been seen in the direction of the Galactic 
center, but the most suitable target for detecting microlensing by Galactic halo objects is the Large 
Magellanic Cloud (LMC), located at galactic longitude and latitude (280°,— 33°) at a distance 
Dgource ~ 50kpc. For source stars in the LMC the characteristic timescale i for microlensing 
events from a lens of mass m is approximately 13Oi/m/M days. The fraction of stars being 
lensed with A > 1.34 is defined as the optical depth, r, which would be [4,10] ~5x 10~ 7 towards 
stars in the LMC if the halo is entirely composed of MACHOs. 

Since July 1992, we have carried out photometric monitoring of about 8.6 million stars in the 
LMC using the 1.27-m telescope at Mount Stromlo Observatory, Australia [11]. 

The data used here consists of our Year 1 data, namely 5169 images of the LMC taken 
between 21 July 1992 and 3 September 1993, with an exposure time of 300 seconds. The images 
are distributed between 22 fields near the center of the LMC [12]. The number of images per field 
ranges from 140 to 350. 

Measurements of the relative brightness of the stars are automatically performed on each 
frame. A good-quality image of each field is first chosen as a 'template', and a list of star positions 
and magnitudes is generated. All other images are aligned with the template both in position 
and flux normalization using a set of bright fiducial stars, and a fit is made to the flux of each 
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star, using a point spread function (PSF) determined from the fiducial stars. Each photometric 
measurement is output with an error estimate and six quality parameters describing the goodness 
of PSF fit, the flux contamination from nearby stars, and the fractions of flux masked by defective 
pixels and cosmic-ray events in the CCD array. 

The dual passband time series for each star are then searched for events consistent with grav- 
itational microlensing. We use empirically determined cuts on the photometry quality flags, along 
with the associated data such as seeing and sky brightness for each image, to exclude questionable 
data points, and we also exclude very red stars which are often irregular variables. The remaining 
lightcurves are then convolved with a set of filters and any lightcurve showing a peak at a modest 
significance level is defined as a 'trigger'. For these triggers, a 5-parameter least-squares fit to 
microlensing is made, defined by Tsit) = A(i)ToB and Fr^) = A{i)TaR, where A(t) is given 
in Eq. 1 above. The fit parameters are the time of peak amplification, t max , the characteristic 
timescale, t, the peak amplification, ^4 max = A(u m \ n ), and the baseline flux of the star in red and 
blue passbands Tqr and Tqb- For each trigger, a large set of associated statistics is also computed, 
and interesting events are automatically selected on the basis of these parameters [13]. 

There are 4 lightcurves from our year 1 data that pass the above cuts: two of these represent 
the same star detected independently in a field overlap, thus we have 3 events consistent with 
microlensing, shown in Fig. 1. The fit parameters are listed in Table 1 below. 



Event 


RA (2000) 


Dec (2000) 


V 


V-R 


*max (days) 


i (days) 


A 


x 2 


1 


05 14 44.5 


-68 48 00 


19.6 


0.6 


57.16 ±0.02 


34.8 ±0.2 


7.20 ± 0.09 


1.34 


2 


05 22 57.0 


-70 33 14 


20.7 


0.4 


121.62 ±0.3 


19.8 ± 1.3 


1.99 ±0.06 


1.41 


3 


05 29 37.4 


-70 06 01 


19.4 


0.3 


154.8 ±0.9 


28.2 ± 1.7 


1.52 ±0.03 


1.01 



Table 1: Parameters of the events. Columns 4 & 5 show approximate magnitude and color 
of the lensed stars. Columns 6-8 show the parameters of the best-fit microlensing models: time of 
peak amplification (Julian days-2449000), the event duration t, and the peak amplification factor, 
with the formal one sigma errors (derived from the covariance matrix of the fit). Column 9 is the 
X 2 per degree of freedom for the microlensing fit. 



Quantitative interpretation of our results depends upon the event detection efficiency S(t) of 
the experiment. Inefficiencies arise because of 1) incomplete sampling of the light curves (primarily 
due to weather interruptions) and 2) because of the blending of two or more stars in an image, 
only one of which will be lensed. In the case of blending, the overall efficiency is reduced because 
of distortion of the light curve (and dilution of the amplification); at the same time the number of 
real target objects is increased, which increases the effective efficiency. These two blending effects 
partially offset each other. 

The efficiency has been determined using a Monte Carlo simulation [13], and the results are 
shown in Figure 2. The upper curve shows the sampling efficiency (neglecting blending effects), 
and is an upper bound. The lower curve includes a treatment of the blending efficiency based upon 
the luminosity function of 'stars' seen in an uncrowded field. The middle curve is based upon a 
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Event 3: A m „= 1.52 t = 28.23 
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Figure 1: The three observed stellar lightcurves that we interpret as gravitational mi- 

crolensing events are each shown in relative flux units (red and blue) vs time 
in days. The solid lines are fits to the theoretical microlensing shape using the 
parameters given in Table 1. 

plausible extrapolation of this luminosity function to fainter stars than we can routinely see in the 
data, and is the 'best estimate' used in this paper. 

In order to illustrate the implications of this result we have compared the number and 
timescales of the observed events with predictions of the commonly used spherical model [1] of 
the density of dark matter in the Galaxy's halo, 



PH(r)=po4r4. ( 2 ) 
+ or 



where r is galactocentric radius, ro = 8.5 kpc is the galactocentric radius of the Sun, a = 5kpc is 
the halo core radius and po = 0.0079 M pc~ 3 = 0.30GeVcm~ 3 is the local dark matter density. 
This model gives [10] a Galactic halo mass of 4.1 x 10 11 M within 50 kpc of the Galactic center, 
and predicts a microlensing rate of T = 1.6 x lO _6 (m/M0) -0 5 events/star/yr. After incorporating 
the detection efficiency, the predicted number of detected events N exp for a dark halo entirely 
composed of compact objects with a unique mass m is shown in Figure 3. Since we have detected 
3 events, models predicting in excess of 7.7 detected events are excluded at 95% confidence. Since 
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Figure 2: MACHO Year 1 microlensing event detection efficiency. Here £ is shown as a 

function of t, the event duration. The left axis is labeled with fractional event 
detection efficiency per monitored object, and the right axis gives the efficiency- 
corrected exposure. The upper curve takes into account only the temporal sam- 
pling effects while the lower two correct for event degradation from blending by 
using two different stellar luminosity functions (LFs) to estimate the increased 
number of targets. The best estimate is the middle curve. 

we are using only the total number of events, note that exclusion of a delta-function MACHO mass 
over a given range also applies to any arbitrary mass function of MACHOs which is contained 
within the same range and provides the same total density. Thus, we see from Fig. 3 that lensing 
objects in the mass range 8 x 10~ 5 M < m < 0.3 M comprise < 100% of the model halo at 95% 
CL. 

Current estimates [14] of the Galactic mass within 50 kpc range from 4 to 6 x 10 11 M , while 
the model of equation 2 contains 4.1 x 10 11 M . Using the 95% CL limits on N exp we can set a 
95% CL upper limit on the mass in this halo of M/ im = 4.1 x 10 11 M@(7.7/N exp ). As shown in 
Figure 3, objects in the mass range 3 x 10~ 4 M < m < 0.06 M contribute no more than 50% of 
the model halo's mass within 50 kpc. 

We have explored a range of different halo density profiles [13], and we find that while the 
constraint on the halo mass fraction in MACHOs is quite model-dependent, our constraints on the 
total mass of MACHOs interior to 50 kpc are relatively independent of the assumed model of the 
galactic halo. 

In setting the limits above we need not assume that our three events are due to microlensing 
by halo objects; however, if we add this assumption we can use the observed durations of the events 
in a maximum likelihood analysis to find the most likely MACHO fraction /, and mass m (again 
with a delta function mass distribution). In this model the rest of the halo is presumed to consist 
of objects with masses outside our range of sensitivity. Figure 4 shows the likelihood contours with 
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Figure 3: The lower panel shows the number of expected events predicted from the stan- 

dard model halo with a delta function mass distribution. Given three observed 
events, points above the line drawn at N exp = 7.7 are excluded at the 95% CL. 
The upper panel shows the 95% CL limit on the halo mass in MACHOs within 
50 kpc of the galactic center for the model. Points above the curve are excluded 
at 95% CL while the line at 4.1 x 10 11 M shows the total mass in this model 
within 50 kpc. 

most likely values of / = 0.20lg 14 an d m = 0.06lg'o4 M© (errors given are the total extent of the 
68% contours). 

The model of equation 2 predicts an optical depth of r TO0( j e i = 4.7 x 10 -7 . Using / we can now 
state the model best fit optical depth of 9.2+g 5 7 x 10 -8 and observed MACHO mass within 50 kpc of 
8.OlJ 4 xlO lo M . We can also estimate the optical depth directly using our observed timescales, our 
exposure E = 9.72 x 10 6 star years, and efficiencies by computing r est = -f^ £(ij/£j) = 8.0 x 10~ 8 , 
in good agreement with the maximum likelihood estimate. 

For three detected events the 90% CL lower bound on the underlying rate is 1.1 events. We 
have estimated [13] the event rates due to the galactic disk, the spheroid and the LMC disk, 
and find that they should contribute on average ~ \ an event to our sample. Our continuing 
observations should help to ascertain the significance of this difference. We caution that recent 
microlensing results toward the Galactic bulge suggest that the standard galactic models used here 
may be incorrect [7-9]. 
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Figure 4: Likelihood contours for the model halo with a delta function mass distribution. 

The numbers labeling the contour levels indicate the total probability of the 
enclosed parameter region according to a Bayesian analysis assuming a prior 
distribution uniform in m and /. The position of the most likely MACHO mass 
m and MACHO halo fraction / is marked with a +. 
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